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Abstract. 

A formalism is presented to evaluate the Sivers function in constituent quark models. A non- 
relativistic reduction of the scheme is performed and applied to the Isgur-Karl model. The sign 
for the u and d flavor contributions that we obtained turns out to be opposite. The Burkardt Sum 
Rule is fulfilled to a large extent. After the estimate of the QCD evolution of the results from the 
momentum scale of the model to the experimental one, a reasonable agreement with the available 
data is obtained. 
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The partonic structure of transversely polarized nucleons is still an open problem[[ll]. 
Semi-inclusive deep inelastic scattering (SIDIS) is one of the proposed processes to 
access the parton distributions (PDs) of transversely polarized hadrons. SIDIS of un- 
polarized electrons off a transversely polarized target shows "single spin asymmetries" 
(SSAs) [120, due to two physical mechanisms, whose contributions can be distinguished 
fiBEl], i.e. the Collins[2] and the Sivers[6] mechanisms. The Sivers mechanism leads 
to a SSA which is the product of the unpolarized fragmentation function with the Sivers 
PD. The latter describes the number density of unpolarized quarks in a transversely po- 
larized target: it is a time-reversal odd. Transverse Momentum Dependent (TMD) PD. 
From the existence of leading-twist Final State Interactions (FSI) 13,181], a non- vanishing 
Sivers function has been explained as generated by the gauge link in the definition of 
TMDs iSH, [O], whose contribution does not vanish in the light-cone gauge, as hap- 
pens for the standard PD functions. Different parameterizations of the available SIDIS 
data have been published [l^, [iH, [13] > still with large error bars. Since a calculation 
from first principles in QCD is not yet possible, several model evaluations have been 
performed, e.g. in a quark-diquark model llv, 15]; in the MIT bag model [16]; in a light- 
cone model [17]. 

We here describe a Constituent Quark Model (CQM) calculation of the Sivers function 



IISQ. CQM calculations of PDs are based on a two steps procedure 1| 191] . First, the matrix 
element of the proper operator is evaluated using the wave functions of the model; then, 
a low momentum scale, /Iq, is ascribed to the model calculation and QCD evolution is 
used to evolve the observable calculated in this low energy scale to the scale of DIS 
experiments. Such procedure has proven successful in describing the gross features 



FIGURE 1. The contributions to the Sivers function in the present approach. 



of PDs and GPDslil, by using different CQMs, e.g. the Isgur-Karl (IK) model 
Besides the fact that it successfully reproduces the low-energy properties of the nucleon, 
the IK model contains the one-gluon-exchange (OGE) mechanism[2l'] . In the present 
calculation, the leading twist contribution to the FSI has to be taken into account. 
We, here, consider the leading, OGE, order, which is natural in the IK model. The 
other approximations in our approach are that, first, only the valence quark sector is 
investigated; second, that the resulting interaction is obtained through a non-relativistic 
(NR) reduction of the relevant operator, according to the philosophy of constituent quark 
models [21]. The Sivers function (Fig. 1) for a proton polarized along the y axis and for 
the quark of flavor =S takes the form 

ft/{xM) = 3!^-ig^^Jdk,dk3^d{k+-xP+)d{k,T+qT-kT)^'^y) 
where g is the strong coupling constant, M the proton mass, and 

^u{d) ^ J- ^l^g^^(k3,m3;ki,mi;P-k3-kum„^ 

^-^-Y^^NR (ki ,k-i,q) 

^sf,s,=- 1 (ki + tn's ■,ki-q,rn[;P-k3~ki, m„ j . (2) 

Using the spin-flavor wave function of the proton in momentum space, ^^f, correspond- 
ing to a given CQM, the Sivers function, Eq. ([T]), can be evaluated. From Eq. Q, one 
notices that the helicity conserving part of the global interaction does not contribute to 
the Sivers function. Besides, in an extreme NR limit, it turns out to be identically zero: 
In our scheme, it is precisely the interference of the lower and upper components in 
the four-spinors of the free quark states which leads to a non-vanishing Sivers function. 
This holds even from the component with / = of the target wave function. While, in 



other approaches 1160, these interference terms arise due to the wave function, they are 
produced here by the interaction. 

The above-described formalism is now applied to the IK model. To evaluate numer- 
ically Eq. ©, g (i.e. (Xs{Q^)) has to be fixed. The prescriptionL19] is used to fix jUq, 
according to the amount of momentum carried by the valence quarks in the model. 
Here, assuming that all the gluons and sea pairs in the proton are produced perturba- 
tively according to NLO evolution equations, in order to have ~ 55% of the momentum 
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FIGURE 2. Left (right): the quantity f^j^^"''''\x), Eq. ©. Dashed curve: IK at Full curve: the 
evolved distribution at NLO. Patterned area: parameterization bv lfisll (see text). 



carried by the valence quarks at a scale of 0.34 GeV^ one finds that /Iq ^ 0. 1 GeV^ 
if Ag^g ~ 0.24 GeV. This yields as{^l)/{An) ~ 0.13 [19]. The results of the present 
approach for the first moments of the Sivers function, defined as 

ft^'^''{x)=jd%^jtf{xM). (3) 

are given by the dashed curves in Fig. 2. They are compared with a parameterization 
of the HERMES data, taken at = 2.5 GeV^ : The patterned area represents the 
1 — a range of the best fit proposed in Ref. [13]. The magnitude of the results is close 
to that of the data, although they have a different shape: the maximum (minimum) is 
predicted at larger values of x. Actually /ig is much lower, = 2.5 GeV^. A proper 
comparison requires QCD evolution of TMDPDs, what is, to large extent, unknown. We 

nevertheless perform a NLO evolution of the model results assuming, for /^^^"^{x), 
the same anomalous dimensions of the unpolarized PDFs. From the final result (full 
curve in Fig. 2), one can see that the agreement with data improves dramatically and the 
trend is reasonably reproduced at least for x > 0.2. Although the performed evolution 
is not exact, the procedure highlights the necessity of evolving the model results to the 
experiment scale and it suggests that the present results could be consistent with data, 
still affected by large errors. 

Properties of the Sivers function can be inferred from general principles. The Burkardt 
Sum Rule (BSR) [|23ll states that, for a proton polarized in the positive y direction, 
I^=.,rf(^f)=Owith 

{kf) = - 1^ dx I dkT^ptf{xM) , (4) 

and must be satisfied at any scale. Within our scheme, at the scale of the model, it 
is found {k") = 10.85 MeV, (fc^) = —11.25 MeV and, in order to have an estimate 
of the quality of the agreement of our results with the sum rule, we define the ratio 



r = I {ki) + {k'^\/\{ki) - {k^) I obtaining r ~ 0.02, so that we can say that our calculation 
fulfills the BSR to a precision of a few percent. One should notice that the agreement 
which is found is better than that found in other model calculations |15, T^, especially 
for what concerns the fulfillment of the Burkardt Sum Rule. However, in a recent work 
[|24ll . we have shown the encouraging result that the calculation in the bag model satisfies 
the Burkardt sum rule at a 5% level. 
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